A dynamical approach of the microcanonical ensemble 
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Abstract 

An analytical method to compute thermodynamic properties of a given Hamiltonian system is proposed. This method 
combines ideas of both dynamical systems and ensemble approaches to thermodynamics, providing de facto a possible alternative 
to traditional Ensemble methods. Thermodynamic properties are extracted from effective motion equations. These equations 
are obtained by introducing a general variational principle applied to an action averaged over a statistical ensemble of paths 
defined on the constant energy surface. The method is applied first to the one dimensional /3-FPU chain and to the two 
dimensional lattice ip 4 model. In both cases the method gives a good insight of some of their statistical and dynamical 
properties. 



The problem raised by Clausius and the second prin- 
ciple found its answer with Boltzmann and the rise of 
equilibrium statistical physics 0, |j. An essential point 
in the theory is related to the law of large numbers, 
which ensures that fluctuations around mean values of 
the thermodynamic quantities are neglig ible H). The 
concept of ensembles is introduced, as for instance the 
microcanonical ensemble for isolated systems, and their 
associated measures are used to average. Developments 
within the Ensemble framework have generalized the use 
of various techniques such as perturbation expansions, 
mean field approximation, or renormalization group Q 
and greatly improved our understanding of phase tran- 
sitions phenomena (see for instance the review H and 
references therein). However, the computation of ther- 
modynamic properties for a given Hamiltonian system 
remains in general inextricable. 




FIG. 1: Representation of the constant energy surface, 
ferent labeled paths are drawn on it. 
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The purpose of this Letter is to introduce an analyt- 
ical approach of the thermodynamic limit and provide 
an alternative to classical techniques. This method relies 
on the large size limit and the universality of trajectories 
(good ergodic and mixing properties are assumed). We 
define an ensemble of paths drawn on the energy surface 
and compute thermodynamic variables through averaged 
equations of motion. This approach applies to systems 
at equilibrium, and proves to be very successful in the 
chosen examples. Note that, the ensemble averaging im- 
plies a large time limit before the large system limit, but 
we invert the order of these two limits. 

Let us identify a set of trajectories on the hypersurface 
defined by the microcanonical measure in the phase-space 
by a set of labels £, which may be initial conditions for in- 
stance. The thermodynamic state does then not depend 
on these labels (this property permits the introduction of 
the ensemble averaging). In the same spirit, we consider 



a family of paths qe(t) (we noted explicitly the time t 
and label I dependences) drawn on the constant energy 
surface (see Fig. |l|). To each path we associate a La- 
grangian L(qi,qi) where the dot denotes time derivative 
and the corresponding action A = J dtL. The basis of 
the proposed method relies on the following claim: since 
the thermodynamic state is label independent, we can 
average the Lagrangian over the labels, and apply the 
variational principle on the mean dynamical system: 



(8A) = 5(A) = 5 / dt (L) = 
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(1) 



(where (• • •) denotes averaging over the labels). The sec- 
ond equality in (|l|) is imposed as a compatibility condi- 
tion at equilibrium and defines a smooth path as the aver- 
age of a flow of paths of the original system. We note that 
after the average is performed, trajectories and points re- 
lated to the mean dynamical system must already com- 
prise some information on the thermodynamic state, 
hence we shall refer to the resulting motion equations 
as thermodynamic motion equations. Let us now con- 
sider Hamiltonian systems of the type H = p 2 /2 + V(q), 
namely quadratic in momentum and with separated con- 
jugated variables. Microcanonical statistics leads to a 
linear relation between the mean kinetic energy MCp 2 /2 
and the temperature T (MC • • • stands for microcanonical 
averaging) m and predicts that the momentum is Gaus- 
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sian with each component pi independent and a vari- 
ance proportional to the temperature MCp 2 ~ T. In 
the canonical ensemble this results in a trivial factoriza- 
tion of the partition function, all the complexity being 
included in the potential V. The present approach uses 
this Gaussian property and reverse the usual argument 
to pass from time averaging to ensemble averaging: at 
thermal equilibrium, we interpret p as being a Gaussian 
stochastic process on the labels and get thermodynamic 
quantities from the mean dynamical system. We now 
propose a possible implementation of these ideas. 

We consider a lattice (in dimension D — 1,2) of N 
sites with coordinates Xi, i = 1, • • • , N. At each site i 
is placed a particle, in general coupled to its neighbors, 
having momentum pi and conjugate coordinate q^. We 
take units such that the lattice spacing, the Boltzmann 
constant, and the mass are equal to one. Since pi is 
Gaussian, we choose to represent it as a superposition of 
random phased waves: 



Nkn 



Pi = E " fe cos(kx l + <j>k) , 



(2) 



k=Q 



where the wavenumber k is in the reciprocal lattice (an 
integer multiple of feo = 27r/A r ( 1 /- D '), the wave amplitude 
is a.},, and its phase 4>k is uniformly distributed on the 
circle. The momentum set is labeled, using (||), with 
the set of phases I = {</>&}. This equation can also be 
interpreted as a change of variables, from p to a, with 
constant Jacobian (the change is linear and we chose an 
equal number of modes and particles). Besides, if the 
total momentum is conserved, we choose to take do = 0. 
As the variance of pi is fixed, we shall assume that the 
dfe are all of the same order (we need a large number 
of relevant modes for the center- limit theorem to apply) . 
Using the relation (p 2 ) = yjd 2 /2 (we average over the 
random phases) and imposing that at equilibrium the 
fluctuations are small, we write (pf) s» T and obtain 
d| rs 0[(T/N)} (we call this relation the Jeans condition 
0). We shall see in the examples that for this scaling 
in N for d and the short range interaction, the mean 
dynamical system becomes a set of oscillators with mean- 
field type interactions and a kind of Jeans spectrum. The 
coordinate variables associated with the representation of 
momenta (||) are 



Nkn 



q. t = a + ^2 a k cos(kXi + <f> k ) 



(3) 
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Note that this equation supposes true the relation pi = 
<ji. The equilibrium state is constructed from the av- 
eraged Lagrangian C = (L)/N, the condition that the 
paths belong to the energy surface e(T) = E/N — 
(H}/N, and the Jeans condition which fixes the tempera- 
ture from the averaged kinetic energy. We in fact applied 
a version of this method to the Kosterlitz-Thouless phase 
transition in the XY model 



We shall start to test this approach with the generic 
case of a chain of coupled harmonic oscillators. The 
Hamiltonian writes H = (1/2 ) V^ -lp 2 + (<7i+i — Qi) 2 ]- 
Using the expressions (||) and (|3|) we compute the av- 
eraged Lagrangian C = (1/4) XLJ"I - u lk a k\ where 
uo 2 )k = 4 sin 2 k/2 and extremize the action to obtain the 
thermodynamic motion equations a k = — w 2 k a k . Equi- 
librium is imposed by the Jeans condition which gives 
a 2 = 2T/Nu 2 . k and leads to the thermodynamic func- 
tion e(T) = T. 

We consider now the Fermi-Pasta-Ulam problem of a 
one-dimensional /3-FPU chain of oscillators. The ther- 
modynamics of this model has been exactly computed 
within the canonical ensemble g]. The Hamiltonian 
reads H = (1/2) J2M + V(q i+ i - g t ) + V{ Qi - q^)}, 
where V{x) — x 2 /2 + [3x A /4. The averaged potential 
energy density v — (V)/N is then 



v = E ( s l a l + 2^ a * ) + 3/3 E hiatal , (4) 

k ^ ~ ' k^k' 

where we used the notation s k = sin(fc/2). It is worth 
noticing that by considering a as the dynamical vari- 
ables of a new system, the interaction v is of the mean- 
field type as the second term in ([|) involves interaction 
between all a k oscillators. The thermodynamic motion 
equations are obtained as before from the Lagrangian 
density £(d, a) = EJ"?/ 2 - v(a)]: 



OLk = 
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Q - 3(3s 2 k a 2 k ] a k , Q 



6/?E4«L (5) 



where Q is a mean-field (intensive) variable. Its fluctu- 
ations at equilibrium are of order 0(l/y/N). Hence we 
consider it in (0) as a constant. This is the large system 
limit N — ► oo taken before the t — > oo limit. In this ap- 
proximation (|5|) describes a set of uncoupled oscillators. 
Moreover, as we verify a posteriori, to the same order of 
approximation we neglect the third term in the brackets 
(it is smaller than the Q term, by a factor 0(1/N)). We 
then obtain a simple linear wave equation with a disper- 
sion relation uj 2 = uj 2 k (l + Q) |1C|| . The Jeans condition 
gives a 2 ujl rj 2T/N, and allows an estimation of the 
neglected term s 2 k a\ rj T/[2N{1 + Q)], Using now the 
dispersion relation, the Jeans condition, and the defini- 
tion of Q, we obtain Q(l + Q)/3/3 = T and the function 
Q = Q(T). Since v = (2Q + Q 2 )/12(3 we finally get the 
thermodynamic relation 



T VI + - 1 

7 + 24/3 



(6) 



This result is compared with the canonical one (||) in 
Fig. |[ The two results are in very good agreement and 
we speculate it is exact. 

For the last example, we consider a two dimensional 
system (D = 2) exhibiting a second order phase transi- 
tion, the so called dynamical lattice <p model studied in 
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FIG. 2: Potential energy versus temperature for = 0.1, the 
dashed line corresponds to the result obtained by jij, while 
the solid line corresponds to equation o. 
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FIG. 3: Solutions of the implicit equation jli]), for m 2 = 
2 and A = 0.6. We notice that depending on the value of 
T up to three different solutions for Q exists giving rise to 
three different branches respectively labelled Qi, Qi, Q3. The 
critical temperature T c is identified. 



Ll|. This model is defined by the Hamiltonian 
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*) 2 , (7) 



where m and A are real parameters, J = 1 is the coupling 
constant, and denotes the summation over the close 
neighbors on a square lattice. In contrast with the other 
cases, where the k = mode ao was a free parameter, in 
this example it is relevant and corresponds to the aver- 
age of q which is proportional to the magnetization of the 
system and is independent of time at equilibrium. The 
computation of the averaged potential energy is very sim- 
ilar to the /3-FPU case: 
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The thermodynamic motion equations are 



= ao ( -a + -Q-m 



(8) 



(9) 
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m + — a n 
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-Q a k + -at , (10) 



where uj 2 k = 4 (sin 2 k x /2 + sin ky/2) is the free harmonic 
frequency (k — (k x , k y ) is now a vector in the plane), and 
Q = a k 1S an intensive variable. 

Equation (^) has multiple solutions in ao depending 
on the temperature through Q — Q(T). Since ao = (qi) 
is the order parameter, we anticipate the existence of a 
phase transition in the thermodynamic state. Indeed, 
the only solution is ao = for Q > = 4to 2 /A but for 



Q < Q* other solutions with finite values of the order 
parameter exist: 



6m 2 3 
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- oQ= o(Q.-Q) 



(11) 



To solve dlQ), we neglect the a\ term, as we did for the (3- 
FPU case (large system limit) and obtain a wave equation 
with the dispersion relation: lo 2 = uo 2 k + fl 2 (T). fl(T) is 
given by, 



O 2 = -(Q.-Q), Q < Q„ 
ft 2 - j(Q-Q*)> Q > Q* 



(12) 
(13) 



where we used the definition of Q* and Eq. (11). We no- 
tice that fl is the frequency corresponding to the absent 
k = mode. 

Given the Jeans condition: a 2 = 2T/ (Nuj 2 ) we notice 
that as long as ft 7^ (Q ^ Q*) the a^ term is a poste- 
riori negligible, we can therefore expect results obtained 
with this approximation to be accurate everywhere but 
at (near) the transition. Another consequence of this 
condition is that most of the modes have comparable 
amplitudes for (1^0, which then implies a Gaussian- 
like distribution for q. At the same order of approxima- 
tion, in the thermodynamic limit, we can compute Q(T), 
identifying it to a Riemann integral. Using the a k given 
by the Jeans spectrum we obtain the following implicit 
equation for Q{T), 



Q(T) = -aK(a) 



n 2 (T) 



(14) 



where K(a) is the complete elliptic integral of the first 

kind rj 2 
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The solutions of Jl4| ) has three branches, plotted on 
Fig. [| We choose the same values m 2 — 2 and A = 0.6 
as the ones used in |Q. We first notice the existence of 
a special temperature close to 19.69 localizes the phase 
transition temperature. Two branches are below Q* = 
40/3, and result from the expression ( |r2j ) used in the 
implicit equation, the third one (on the top of the figure) 
corresponds to the expression (|13f ) . The fact that the 
wave form solutions are not valid for Q — Q* translate 
in the divergence of K(a), since as Q — * Q*, a — ► 1. The 
divergence is logarithmic, and then for sufficiently small 
T, a solution around Q» of jl^ ) always exists; resulting in 
the two upper branches being asymptotic to the Q{T) = 
Q* curve as T goes to 0. 

In order to select one branch from another we compute 
their respective density of energy. According to equations 
and (H), we have two different expressions for the 
density of energy, 

l r -|Q 2 > Q>Qc ' [ ' 

The results are plotted as the temperature versus the 
density of energy in Fig. || in analogy to the results 
presented in [^lj. The physical relevant solution is the 
one whose energy is the smaller for a given temperature, 
which translate in the upper line in the figure. The tran- 
sition is then identified at a density of energy e c = 25.07, 
whose corresponding temperature is T c = 19.69. These 
results are in good agreement with the one predicted by 
numerical simulations |p4| : respectively T c = 17.65 and 
e c = 21.1. Using equation Q), we have also access to 
the square of the magnetization and there is also a good 
quantitative agreement with the numerical results. The 
discontinuous behavior of the magnetization at the tran- 
sition is although surprising. This behavior was also ob- 
served numerically in jn], and is explained by noticing 
that the true order parameter is (\q\) and not \(q}\. In 
the present case, we can also wonder whether this behav- 
ior is due to the N — ► oo limit taken before the t — ► oo, 
Indeed the neglected terms are relevant at the transi- 
tion, and only become negligible around the transition 
after the N — > oo limit, which may affect the nature of 
the observed transition. However this behavior may also 
find its origin in the choice of writing the momentum as 
a superposition of N' random phased waves (motivated 
by the solutions of the linearized equations of motion) 
equal to the number of degrees of freedom N. Writing 
the momentum with the number of modes N' being a 



growing unbounded function of TV is sufficient to obtain 
a Gaussian process. Another representation may then 
be appropriate to tackle the transition region and for in- 
stance in H a high temperature approach was used to 
compute the critical temperature. 

To conclude we point out that the thermodynamic mo- 
tion equations method allowed us to compute the macro- 
scopic properties of coupled nonlinear oscillator systems 
in one and two dimensions. Quantitative agreement with 
exact or numerical results of these quantities is obtained. 
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FIG. 4: Temperature versus density of energy for m 2 = 2 and 
A = 0.6. The three different branches resulting from the solu- 
tions of Q presented in Fig. ^ are represented. The physical 
line is the upper one. The transition is clearly identified with 
T c w 19.64 and e c ~ 25.07. 



Moreover the phase transition for the ip 4 model is de- 
tected and a good estimate of the critical energy and 
temperature are given even though we approximatively 
solved the thermodynamic motion equations. We expect 
that this method will be successful for other systems and 
speculate that the actual solving of the exact thermo- 
dynamics motion equation should lead to an exact ther- 
modynamic limit. We believe it may also be possible to 
extend the scope of the method to systems out of equilib- 
rium and describe their macroscopic evolution with the 
thermodynamic motion equations. 
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